|
In mathematics, the power iteration is an eigenvalue algorithm: given a matrix ''A'', the algorithm will produce a number λ (the eigenvalue) and a nonzero vector ''v'' (the eigenvector), such that ''Av'' = λ''v''. The algorithm is also known as the Von Mises iteration.〔Richard von Mises and H. Pollaczek-Geiringer, ''Praktische Verfahren der Gleichungsauflösung'', ZAMM - Zeitschrift für Angewandte Mathematik und Mechanik 9, 152-164 (1929).〕 The power iteration is a very simple algorithm. It does not compute a matrix decomposition, and hence it can be used when ''A'' is a very large sparse matrix. However, it will find only one eigenvalue (the one with the greatest absolute value) and it may converge only slowly. ==The method== The power iteration algorithm starts with a vector ''b''0, which may be an approximation to the dominant eigenvector or a random vector. The method is described by the recurrence relation : So, at every iteration, the vector ''b''''k'' is multiplied by the matrix ''A'' and normalized. If we assume ''A'' has an eigenvalue that is strictly greater in magnitude than its other eigenvalues and the starting vector has a nonzero component in the direction of an eigenvector associated with the dominant eigenvalue, then a subsequence converges to an eigenvector associated with the dominant eigenvalue. Without the two assumptions above, the sequence does not necessarily converge. In this sequence, : , where is an eigenvector associated with the dominant eigenvalue, and . The presence of the term does not converge unless defined by : converges to the dominant eigenvalue. This can be run as a simulation program with the following simple algorithm: The value of ''norm'' converges to the absolute value of the dominant eigenvalue, and the vector ''b'' to an associated eigenvector. ''Note:'' The above code assumes real A,b. To handle complex; A()() becomes conj(A()()), and tmp() *tmp() becomes conj(tmp()) *tmp() This algorithm is the one used to calculate such things as the Google PageRank. The method can also be used to calculate the spectral radius (the largest eigenvalue) of a matrix by computing the Rayleigh quotient : 抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)』 ■ウィキペディアで「Power iteration」の詳細全文を読む スポンサード リンク
|